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QCD Critical Region and Quark Gluon Plasma from an Imaginary fiB 
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We discuss the imaginary chemical potential approach to the study of QCD at nonzero temperature and density, 
present results for the four flavor model in the different phases and show that this method is ideally suited for a 
comparison between lattice data and phenomenological models. 



1. QCD and a complex /is 

Results from simulations with an imaginary 
chemical potential can be analytically continued 
to a real chemical potential, thus circumventing 
the sign problem P, |2| [3]. In practice, the an- 
alytical continuation is carried out along one line 
in the complex ^ plane: first along the imaginary 
axes, and then along the real one. It is then mean- 
ingful to map this path in the complex /i^ plane: 
because of the symmetry property Z{fi) = Z{—^) 
this can be achieved without losing generality. In 
the complex /x^ plane the partition function is real 
for real values of the external parameter /z^ , com- 
plex otherwise: the situation resembles that of 
ordinary statistical models in an external field. 
Hence, the analyticity of the physical observables 
^ as well as that of the critical line ^ follows 
naturally. 

The phase diagram in the temperature, (real) 
jjL^ plane is sketched in Fig.l, where we omit 
the superconducting and the color flavor locked 
phase, which (unfortunately) play no role in our 
discussion. The region accessible to numerical 
simulations is the one with /i^ < 0: at a vari- 
ance with other approaches to finite density QCD 
which only use information at = [S] [Hj |7| 
the imaginary chemical potential method exploits 
the entire halfspace. And we will also argue that 
there are physical questions which can be ad- 
dressed without analytic continuation. 



Figure 1 . Sketch of the phase diagram in the /i^ , T 
plane: the solid line is the chiral transition, the 
dashed line is the Roberge Weiss transition. Sim- 
ulations can be carried out at < and results 
continued to the physical domain /j.^ > 



2. The Critical Line 

Also in this case the consideration of the T, 
plane helps the analysis. Model analysis suggests 
the following parametrization, confirmed by nu- 
merical results: 



(1) 
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It encodes reality for real /i^, contains the phys- 
ical scale Tc , is dimensionally consistent, gives 
r(/i = 0) = Tc, T(/x 9^ 0) < Tc . We refer to Sec- 
tion IV of Ref. P] for our results on the critical 
line in the four fiavor model, and their discussion 
in terms of model calculations. Here it suffices to 
remind ourselves that the second order approx- 



2 



S 


c o 52 








a 

■ 


• T = 0.95 T 
T = 0.985 
T T = 1.04 < T„ 








- T = 1.065 T, < T, 




8 




■ T = 1.095 = T, 
□ T = 1.5 
i T = 2.5 








T = 3.5T, 








■ 









0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



■ 5 = 



5.07 5.08 5.09 5.1 5.11 5.12 5.13 



Figure 2. Number density as a function of /i/: 
note the smooth behavior in the hadronic phase, 
consistent with the hadron resonance gas model, 
the chiral transition in the Roberge Weiss region, 
the rapid increase in the plasma phase, approach- 
ing a nearly free quark behavior. 



imation turns out to be adequate, and that the 
fourth order corrections were found to be consis- 
tent with zero within errors. 

3. Hadronic Phase: T < Tc 

In this region observables are a continuous and 
periodic function of /z/ /T , analytic continuation 
in the ij? > half plane is always possible, but 
interesting only when = 0, T) > 0. 

The analytic continuation of an observable O is 
valid till fi < ^ic{T), where ^iciT) has to be mea- 
sured independently. The value of the analytic 
continuation at ^Lc, 0{fic) defines the discontinu- 
ity at the critical point, or, equivalently, the crit- 
ical value of, say, the number density. In turns, 
this allows the identification of the order of the 
phase transition. 

For observables which are even/odd in the 
chemical potential, Og/o, we have considered a 
Fourier series, observing that the first cosine/sine 
terms suffice to parametrize the data '3' . This has 
been confirmed in where this result has been 
intepreted within the framework of the hadron 



Figure 3. Correlation between < ipip > and 
Polyakov loop at = 0.15, demonstrating the 
correlation of chiral and deconfining transition at 
nonzero baryon density. 



resonance gas model. 

4. Roberge Weiss Regime: Tc <T <Te 

The analytic continuation is valid till ^ — oo 
but the interval accessible to the simulations at 
imaginary /i is small, as simulations in this area 
hits the chiral critical line for /i^ < 0. 

The bright side of this is that the nature of the 
critical line can then be studied without need for 
analytic continuation. In Fig. 3 we show the clear 
correlation between the Polyakov Loop and the 
chiral condensate at = 0.15. The correlation 
between chiral and deconfining transition persists 
at nonzero chemical potential. 

It is also of interest to note that the non- 
applicability of perturbation theory in this region 
is almost a theorem: indeed the analytic continu- 
ation of the polynomial predicted by perturbation 
theory for positive fi^ would never reproduce the 
correct critical behavior at the second order phase 
transition for /i^ < , and it is then ruled out. 

5. The QGP phase : Te <T 

Several analytic models have been proposed to 
describe the properties of this phase and, as ana- 
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Figure 4. The effective prefactor plot (see text). 



lytic models can be obviously analytically contin- 
ued in the /i^ < half plane, imaginary chem- 
ical potential is an excellent testbed for these 
approaches. Here, as an example, we just con- 
trast the data with a free field behavior AP/T^ — 
2{p./Tf , where AP = P{y) - P(0), and we have 
ignored the fourth order terms. In Fig. 4 we plot 
APK^N^^4,)/{T'^2{fi/T)^) versus fi/T, where we 
have corrected for finite lattice effects Ki^(^]y^^4) 
following 0, El- 

This "effective prefactor plot" is perhaps more 
informative than the direct quadratic fits of 
AP/T"* to fc(/i/r)^, as it allows an assessment by 
eye of the /i dependence, if any, of the prefactor 
to the quadratic term. 

We see that the results approach the perturba- 
tive limit, but corrections at small chemical po- 
tential are clearly visible, possibly consistent with 
the predictions of |10| . 

6. Summary/Outlook 

We have studied four flavor QCD within the 
imaginary chemical potential approach in a large 
part of the phase diagram. 

We have shown that the chiral and deconfining 
transition remain correlated at nonzero chemical 
potential. We have found that the critical line is 
described by a polynomial, and interpreted this 



result in terms of simple models. We have iden- 
tified and discussed three different regimes: the 
hadronic phase results are consistent with the 
hadron resonance gas model. The Roberge Weiss 
regime is eminently nonperturbative, and in this 
regime we have the possibility to study the na- 
ture of the chiral transition at nonzero chemical 
potential without performing any analytic con- 
tinuation. The Quark Gluon Plasma phase is an 
ideal test bed for analytic models. 

The method seems mature for quantitative 
studies in realistic cases, and a nice possibility 
is offered by a combination of this approach with 
other methods, for instance by using reweighting 
|H] [7] or direct calculations of derivatives at 
nonzero fi to improve the accuracy of the results 
at negative /i^ . Finally, the study of discontinu- 
ities as sketched in Sect. 3 above might offer an al- 
ternative approach to the study of the endpoints 
and tricritical points. 
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